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■ We have obtained optical spectra of 29 early-type (E/SO) galaxies that hosted 
j_] ■ type la supernovae (SNe la). We have measured absorption-line strengths and 
^ compared them to a grid of models to extract the relations between the super- 
nova properties and the luminosity-weighted age/composition of the host galax- 
ies. Such a direct measurement is a marked improvement over existing analyses 
which tend to rely on general correlations between the properties of stellar pop- 
ulations and morphology. Our galaxy sample ranges over a factor of ten in iron 
abundance and shows both old and young dominant population ages. The same 
analysis was applied to a large number of early-type field galaxies selected from 
the SDSS spectroscopic survey. We find no difference in the age and abundance 
distributions between the field galaxies and the SN la host galaxies. We do find 
a strong correlation suggesting that SNe la in galaxies whose populations have 
a characteristic age greater than 5 Gyr are ~ 1 mag fainter at Vmax than those 
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found in galaxies with younger populations. However, the data cannot discrimi- 
nate between a smooth relation connecting age and supernova luminosity or two 
populations of SN la progenitors. 

We find that SN la distance residuals in the Hubble diagram are corre- 
lated with host-galaxy me t al ab undance, consistent with the predictions of 



Timmes. Brown, fc Truraru (120031 ). The data show that high iron abundance 
galaxies host less-luminous supernovae. We thus conclude that the time since 
progenitor formation primarily determines the radioactive Ni production while 
progenitor metal abundance has a weaker infiuence on peak luminosity, but one 
not fully corrected by light-curve shape and color fitters. This result, particularly 
the secondary dependence on metallicity, has significant implications for the de- 
termination of the equation-of-state parameter, w = P/^pc^), and could impact 
planning for future dark-energy missions such as JDEM. 

Assuming no selection effects in discovering SNe la in local early-type galaxies, 
we find a higher specific SN la rate in E/SO galaxies with ages below 3 Gyr than 
in older hosts. The higher rate and brighter luminosities seen in the youngest 
E/SO hosts may be a result of recent star formation and represents a tail of the 
"prompt" SN la progenitors. 

Subject headings: galaxies: supernovae and spectroscopy - cosmology: distance 
scale - supernovae: general 



1. Introduction 



Type la supernovae (SNe la) have proven to be effective standardizable candles and their 
high peak luminosity has made them excellent cosmological probes. An em pirically deter- 



mined relation between light-curve properties and peak absolute magnitude (IPhillips 



Hamuv et al.lll996bl: iRiess. Press. &: Kirshnerlll996l : iGoldhaber et al.l l200ll : iGuy et al 



Jha. Riess. &: Kirshnei 



Hubble constant (jjha et al. 



200 



3 



Conlev et al 



1999 



1993'; 
2007i : 



2008h has yielded improved estima t es of the global 



Freedman et al.l l200ll : iRiess et al.l l2005l : ISandage et al. 



20061 ) . Although the earliest atte mpts to measure the matter de nsity of the universe preceded 
these light-curve shape methods (INorgaard-Neilsen et al.lll989l ). and the initial efforts by the 
Superno va Cosmology Project ( SCP) led to incorrect conclusions about the cosmic matter 
density (IPerlmutter et al.l 119971). earlv work by both the SCP and the High-Z Supernova 



Search Team (ISchmidt et al. 



1998: Garnavich et al 



1998aJ) 



19981 ) found evidence for low matter density (IPerlmutter et al. 



Subsequent publications by iRiess et al.l (119981 ) and 



Perlmutter et al 



([19991) came to 
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the surprising conclusion that the universe is exp a nding at a n accele ratin g rate, driven 



by a mysterious dark energy; see iKirshnerl (|2002| ). iFilippenkol (l2005bl ). or iFrieman et al. 
( 120081 ) for re views. This surp r ising result has been confirmed by more-recent super nova 
observations (iKnop et al.ll2003l : iTonry et al.ll2003l : iBarris et al.ll2004l : iRiess et al.ll2004j ) and 



by co mplementary measu rements of the cosmic microwave background (CMB) anisotropies 



Bennett et al.ll2003l ) and large-scale structure (e.g.. iPeacock et al.ll200ll ). Focus has 



(e.g. 

now shifted from demonstrating the existence of dark energy to constraining its proper 



ties ( 


Garnavich et a . 


1998bl: 


Astier et al. 


200d: 


Riess et a . 


2007: 


Wood-Vasev et al. 


2007; 


Miknaitis et al. 


2007; 


Eisenstein et al.l 


2OO7I; 


Frieman et al. 


20071) and sharDenine the 


tools 


available for its study (e.g., 


Krisciunas et al. 


2007; 


Wood-Vasev et al. 


2008). 



Analyses of their spectra and light curves have led to the consensus that SNe la are well 
characterized by a thermonuclear disruption of a carbon-oxygen white dwarf (WD). The 
transition from a dynamically stable WD to an explosion capable of outshining an entire 
galaxy is precipitated by the WD mass approaching the Chandrasekhar mass limit through 
matter accretion from a close binary companion. The nature of this accretion process is 
still uncertain, though two scenarios have become the focus of SN la progenitor model- 
ing. The first is the single- degenerate (SD) mode l in which a W D accretes matter from a 
non-degenerate companion (jWhelan fc Ibenl Il973l : iNomotd Il982l ). while the second is the 



double-degener ate (DD) model involving the merger to two WDs (llben fc Tutukovl Il984 



Webbinklll984l ). Though promising, our understanding of these models remains incomplete 
with questions persisting about the nature of the binary companion (main sequence, asymp- 
totic giant branch, etc.) and accretion flow in the SD scenario, the apparent likelihood 
that a WD merger will result in an accretion-induced collapse rather than a SN la, and the 
potential for multiple stellar-population-dependent channels leading to SN la explosions. 

Differentiating among dark-energy models requires mapping the expansion history of the 
universe with high precision and a small systematic error. Applying SNe la to this problem 

f the explosion physics and progenitor stars (e.g.. 



Howell et al. 


2006; 


Hicken et al. 


2007) 



results from the thermonuclear disruption of a near- Chandrasekhar mass carbon-oxygen 
WD, there remain p otential dangers iri blindly applying SNe la to the dark-energy problera . 



Theoretical models (IHoflich et al.lll998l : lUmeda et al.lll999l : iTimmes. Brown, fc Truranll2003l ) 



suggest that population age and metallicity may have an effect on the WD composition and 
influence the SN la peak luminosity. Since the average age and metal content of stars have 
evolved over cosmic time, the ensemble character of SN la explosions may be a function of 
lookback time. But the local universe contains a range of stellar ages and metallicities, so 
an empirical calibration of these effects may be possible. 
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Hamuy et al.l (Il996bl ) were the first to note a correlation between SN la decline rate 
and host-galaxy morphology. The results showed that intrinsically faint events occur in 
early-type (E/SO) galaxies, while luminous events are often hosted by late-type galaxies. 
This correlation with morphology has since been confirmed in larger SN la samples (e.g., 
Gallagher et al.ll2005l ). but the cause of the trend remains unknown. The existence of a 
correlation is important because it demonstrates that the SN la environment imprints itself 
on some aspect of the progenitor. On average, spiral and elliptical galaxies differ in popu- 
lation age, star- formation rate (SFR), and metal conten t , with significant overlap in these 
properties across the morphological types. iHamuy et al.l (|2000| ) measured the metallicity of 
a small sample of star-forming hosts to test if age or metal content was the factor influencing 
SN la luminosity, but the results were ambiguous. 



Gallagher et al.l (120051 ) expanded the number of observed hosts to test the impact of 
stellar environment on the photometric properties of nearby SNe la. Potential dependencies 
of SN decline rate on host-galaxy absolute magnitude and star-formation history were inves- 
tigated in a sample of galaxies with a broad range of properties and across the full Hubble 
sequence. Metallicity was measured from emission lines, limiting this aspect of the study 
to star-forming hosts. No significant correlation between metallicity and decline rate was 
detected, though a tenuous correlation was found between the metal content and Hubble 
residuals. Furthermore, only galaxies without significant star formation were found to host 
faint, fast-declining SNe la. The division by star-formation history was more stark than by 
morphology alone, suggesting that the observed luminosity spread in SNe la results from the 
influence of population age on the WD progenitors. The luminosity /star-format ion relation 
has also been s een in a high-redshi ft sample of SNe la studied by the Supernova Legacy 
Survey (SNLS) JSuUivan et al.ll2006h . 



The specific SN la rate is strongly dep endent on the amount of current star formation 
( iMannucci et al.l l2005l : ISuUivan et al.l l2006l ) . Actively star- forming galaxies host approxi- 
mately 10 times as many SNe la per unit mass than do their low star- forming counterparts. 
Scannapieco fc BildstenI (120051 ) find that the SN la rate is well described by a combination 
of two distinct populations of SNe la: a "prompt" (shorter delay time) SN la component 
dependent upon the recent star formation in the galaxy and an "delayed" or "tardy" (longer 
delay time) component dependent on the number of low-mass stars. Assuming the light- 
curve shape correction methods are perfect, then this will not have an effect on the recent 
cosmological results. However, if the corrected photometric properties of these two SN pop- 
ulations differ (e.g., a variation due to age or metallicity), then it could have implications for 
derived cosmological parameters because the short delay time SNe are expected to dominate 
at high redshift whereas the long delay time events are expected to dominate at low redshift. 
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Here we present our study on the luminosity- weighted ages and metal abundances of 
early-type host galaxies. Early-type galaxies, unlike star-forming spirals, have been shown 
to host both faint and bright SNe la, thereby allowing us to sample the full diversity 
of SNe la. Consequently, we focus our study on absorption-line spectra of host galax- 
ies and utilize sing le age, single rnetall i city stellar population models to characterize the 
stellar populations (Vazdekis et al.lll996l : IVazdekislll999l : iBlakeslee. Vazdekis. fc Aihaiil200ll : 
Vazdekis et al. 2003120081 ). The impact of both age and metallicity on SN peak magnitude is 
studied and the sample age and metallicity distributions will be compared to absorption-line 
spectra of field galaxies obtained from the Sloan Digital Sky Survey (SDSS). Furthermore, 
we will infer from our data a measurement of the relative SN la rate as a function of galaxy 
age. Finally, we will search for second-order dependencies at the level of the SN la intrinsic 
scatter by comparing host-galaxy ages and metallicities to the SN residuals from the Hubble 
diagram. 

We detail our observational techniques and data-reduction pipeline in §2 of this paper. 
In §3 we report the methods for performing emission-line corrections, and for obtaining 
age and metallicity estimates of our host sample. Our results are presented in §4, and the 
summary and conclusions are in §5. 



2. Observations and Data Reduction 



2.1. Observations 



The objective of this study was to obtain absorption-line spectra of SN la host galax- 
ies with the intent of estimating their ages and metallicities throug h a comparison with 



the simple stellar population (SSP) models of Alexandre Vazdekis (IVazdekis et al.l 12008 



Sanchez-Blazquez et al.ll2006l ). To this end, we obtained spectra of a sample of SN la host 
galaxies with the 1.5 m Tillinghast telescope located at t he F. L. Whipple Observatory. We 



employed the FAST spectrograph flFabricant et al.lll998l ) fitted with the 300 line mm~^ re- 
flection grating and a 3" x 3' slit. The setup gave a resolution of ~ 7 A full width at half 
maximum intensity (FWHM). 

The spectra were obtained over the course of 6 nights during April and July of 2005. 
Our s et was comp iled f rom the host galaxies of SN e la from the samples of iPhillips et al. 
([19991), LJhal (12002), and iGaneshalingam et al.l (120081 ). Suitable targets were those classified 



as having early type (SO-E), as well as spiral galaxies contain ing strong absorption hues 
in their spectra (and minimal emission lines) as determined by iGallagher et al.l (120051 ). In 
addition to the host galaxies, we also obtained spectra for a set of comparison elliptical 
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galaxies having age and metallicities estimates determined by lTrager et al.l (120001 ) . 



The sht was ahgned along the galaxy's major axis with position angles (PAs) determined 
via the Digital Sky Survey (DSS) plates. In a few cases the PA was altered to prevent a nearby 
star from falling on the slit. Exposure times were varied depending on the galaxy brightness, 
with a target signal-to-noise ratio (S/N) of 10-20. For each target we obtained three spectra 
that we subsequently combined to improve the S/N and remove cosmic rays. Bias, dark, and 
dome-fiat exposures were taken at the beginning and the end of each night, and a comparison 
spectrum of a He-Ne-Ar lamp was acquired before each target for wavelength calibration. 
Finally, flux stan dard star spectra were taken each night with the slit oriented along the 
parallactic angle (lFilippenkdll982l ). 



Our full sample of galaxies can be seen in Table 1. Columns (1) and (2) give the 
host-galaxy name and hosted SN, respectively. Column (3) provides the peak magnitudes 
of the hosted SN la, determ ined from the SN light-curve data and employing MLCS2k2 
( jJha. Riess. fc Kirshnerl 120071 ). Column (4) gives the morphological classification of each 
galaxy, while Column (5) reports the PA for the each observation. Column (6) shows the 
width of the extracted aperture encompassing all the light out to approximately the effective 

gives the source of the SN light-curve: (1) the CfA sam- 



radi us, Rf l|. Finally, Colunin (7 



( Ganeshalingam et a 



( iTaubenberger et al 



2008 



pie (IJha. Riess. fc Kirshnerl 120071), (2) the Katzman Automatic Imagi ng Telescope (KAIT) 



Filippenko et al 



20081 : ICallagher et al.lboOsX 



gOOll : lFilippenkoll2005af ). and (3) SN 2005bl 



2.2. Data Reduction 

The majority of data reduction for this study was performed using standard techniques 
within IRAm The data were bias and dark-subtracted, and flat-fielded to remove the pixel- 
to-pixel variations in the CCD; also, the bad rows and columns of the CCD were interpolated 
across using the routine FIXPIX. Following preprocessing, the task was to combine the set 
of three spectral images that we had for each target. The first step was the removal of cosmic 
rays. In the interest of preserving every image in a set, thereby increasing our statistics, we 



^The radius within which half of the total galactic light is contained. 

^IRAF is the Image Reduction and Analysis Facility, a general purpose software system for the reduction 
and analysis of astronomical data. IRAF is written and supported by the IRAF programming group at the 
National Optical Astronomy Observatories (NOAO) in Tucson, AZ. NOAO is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc. under cooperative agreement with the National 
Science Foundation (NSF). 
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removed each cosmic ray individually using the IMEDIT routine in IRAF. Cosmic rays that 
fell on the galaxy nucleus itself were left alone until after the aperture extraction. 

One-dimensional aperture extraction was performed using the APALL routine with 
extraction radii equal to the effective radius (Re) of each galaxy. The effective radii were de- 
termined through an IDL code written to integrate over the average profile of 400 centrally 
located columns along the spatial axis of each image. Each spectrum was then individu- 
ally checked for cosmic rays landing on one of our relevant absorption-line indices or the 
surrounding continua. In the event that an interfering cosmic ray was found, steps were 
taken to either remove it interactively within SPLOT, or to remove the image altogether 
if the cosmic ray was too disruptive. The average of each set of images was generated us- 
ing IM-COMBINE, and the resultant spectra were wavelength and flux calibrated. Finally, 
the spectra were corrected for Ga lactic extinction using the IRAF routine DEREDDEN 



fICardelli. Clavton. fc Mathislll989[ ). 



3. Disentangling Age and Metallicity 
3.1. Background 

A well known and much reviled puzzle in the study of stellar populations is the observed 
similarity between the effe cts of age and metallicity on the integrated light of stellar popula- 



tions (e.g.. IWortheylll994l ). Astronomers attempt to unravel this age/metallicity degeneracy 
through the development of stellar population models. The simplified goal of stellar popula- 
tion synthesis modeling is to find a combination of stars for which the integrated spectrum 
of the stars matches the observed spectrum of the population under study. 



Early empirical techniques such as quadratic programming devised by lFabeii (119721 ) have 



1970s n 


rinslev & Gunn 1976: Tinslev 1980 


Gunn. Strvker, & Tinslev 1981: Worthev 


1994; 


Buzzoni 


1995: Bressan. Chiosi. & Tantalo 


1996: Maraston 


1998: Bruzual 


& Chariot 


2003,; 


Fioc & Rocca-Volmeraned 1997: Vazauez & Leitherei 2005: 


Vazdekis et al. 


2003: Dekado et al. 


2005). Whereas empirical models are constructed from a combination of stellar spectra or 



the spectra of stellar clusters, evolutionary synthesis models utilize theoretical stellar evolu- 
tionary isochrones as the primary constituent. A single isochrone on the Hertzsprung- Russell 
(HR) diagram gives the locus of luminosities and temperatures at a single moment in time 
for stars of all masses. Together with an assumed initial mass function (IMF), an isochrone 
can be built to accurately model a single-age stellar population. The final model is obtained 
by converting isochrone parameters to observed spectra and finally integrating along the 
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isochrone. 

The populations that are modeled using the above techniques are known as simple stellar 
populations, or SSPs. SSPs are, by definition, a population of stars created during a single 
burst event (i.e., all the stars have the same age) and possessing a single global metallicity. 
However, early-type galaxies are often comprised of multiple-component stellar populations. 
Galaxy ages derived through a comparison with these single-component SSPs are luminosity 
weighted, and are typically sensitive to the youngest component in a real early-type galaxy. 
Consequen tly, it is appro priate to view the ages derived in studies such as ours as lower-limit 



estimates (jRenzinill2006l ) 



The determination of age and metallicity via galactic spectra is fairly simple, in principle. 
We can exploit predictions made by population synthesis models to relate physical properties, 
such as age and metallicity, to observables, such as the absorption-line strengths for the 
strongest atomic and molecular absorption features in the optical range. Measuring the 
relative strengths of these absorption features, we can estimate the age and metallicity of 
the stellar population that created the spectrum. 



3.2. Model Broadening 



For the age and metallicity estimates of our g alaxies, we emp l oy th e single age , singl e 
metalhcity stellar population mode l s described bv IVazdekis et al.l (1996). IVazdekisi (19991 ). 
Blakeslee. Vazdekis. fc A^h^ J200lh . IVazdekis et all J2003h . and lVazdekis et all J2008h. The 



models make use of the e mpirical stellar spectral library MILES (jSanchez-Blazquez et al. 



20061 : ICenarro et al.l 120071 ). MILES is a marked improvement over previous libraries; it has 
985 stars, a spectral resolution of 2.3 A (FWHM), and a wavelength range of 3525-7500 A. 

One advantage of empirical libraries over theoretical stellar libraries is that they are 
based on actual stellar spectra and are not dependent on the potential inaccuracies and 
underlying assumptions inherent in any theoretical model. However, in order for an empirical 
stellar library to sufficiently cover the parameter space in temperature, gravity, abundance, 
and [a/Fe], the stellar observations must be of the highest quality. Consequently, empirical 
stellar libraries have traditionally been restricted to the nearest stars, thereby introducing 
a bias of atmospheric parameters toward those which are seen in the solar neighborhood. 
The MILES library acquired atmospheric parame ters from the literature and calib rated the 
set throug h a subset of r e ferenc e field stars from ISoubiran. Katz. fc Cayrell (119981 ): for full 

(120071 ). In this way, MILES has optimized its stellar atmosphere 



details see Cenarro et al 



coverage and improved upon previous libraries. Nevertheless, because early-type galaxies 
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are generally a enhanced and metal rich, the potential bias toward solar abundance ratios 
should be kept in mind. 

The stellar population synthesis models derived from MILES consist of single age, sin- 
gle metallicity spectral energy distributions (SEDs) that adopt a standard Salpeter IMF 
(ISalpeterlll955l ). The SEDs possess a FWHM resolution of 2.3 A and a spectral range simi- 
lar to our data at 3540-7410 A. There are 276 models with ages ranging from 0.10 to 17.78 
Gyr and metallicities ranging from [M/H] = —1.68 to 0.20. 

In the past, the most widely used method of estimating stellar population age and metal- 
licity ha s been through a comparison of absorption- line indices with those of the Lick/IDS 
System (jWortheyi 1 19941 ) . Since the stars of the Lick/IDS stellar library are not flux cali- 
brated, data must be conver ted to the instrumental response curve of the original data set 
( iWorthey &: Ottavianil 119971 ). With the advent of improved, fiux-calibrated stel lar libraries , 
it became possible to genera t e complete SEDs of a sing le stellar population (jjonesl 119971 : 
Sanchez-Blazquez et al.l l2006l : IVazdekis fc Arimotd 119991 ) , thereby eliminating the need for 
such a correction. Furthermore, in order to use the Lick population models, it is necessary 
to transform the observational data to match the resolution of the Lick spectrograph (~ 8- 
10 A FWHM). For our analysis, the model SEDs are at a higher resolution than the data; 
consequently, the models were required to be broadened to the resolution of our data rather 
than the other way around. 

Thus, the first step in the analysis was to characterize the wavelength-dependent reso- 
lution of our spectra. This was accomplished by measuring the line widths of emission lines 
in the He-Ne-Ar lamp across our full spectral range. A polynomial fit to the resulting scat- 
ter plot provided a wavelength-dependent resolution function for the FAST spectrograph. 
We derived a wavelength-dependent broadening term with a square equal to the difference 
between the square of the host-galaxy resolution and that of the resolution of the Vazdekis 
models, given as 2.3 A. 

The next broadening term we considered was from the velocity dispersion of the host 
galaxies. To determine the velocity dispersion for each galaxy, we obtained spectra of ve- 
locity standard stars HD12623 and HD52071; these were taken within a few months of our 
observing runs on the Tillinghast telescope using the FAST spectrograph. In total, we ob- 
tained 6 spectra of HD12623 and 14 spectra of HD52071. Cross-correlation analysis was then 
performed between these templates and the host-galaxy spectra using the IRAF routine FX- 
COR. A velocity dispersion and redshift estimate were obtained using each template-galaxy 
combination, with the final result for a given host being the average of the measurements 
made with each of the twenty templates. 
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Our two broadening terms were summed in quadrature, yielding the final width of our 
Gaussian kernel used in the convolution with the host galaxy. Each host galaxy had a 
unique final broadening term and the Vazdekis model SEDs were broadened independently 
according to this term. In this way, we created separate age/met allicity grids for each host 
galaxy. The velocity dispersions and heliocentric velocities resulting from this analysis are 
shown in Table [2l These heliocentric velocities were converted to CMB rest frame for the 
subsequent analysis using the NED velocity calculator. 

The age dependent absorption index chosen for this study was H/3. H/5 is known to be 
metallicity insensitive and older stellar populat ions are knoAvn to have relatively strong B.0 



absor ption, eliminating the need for high S/N (jCardielll2003l : ICaldwell. Rose, fc Concannon 



20031 ). Although our sample was comprised mostly of older stellar populations, there was a 
chance for H/3 emission contamination. Consequently, our data were emission-line corrected 
with the procedure discussed in the next section. We chose two age- insensitive Fe absorption 
lines for our abundance analysis, FeA5270 and FeA4383. Plotting both of these iron lines 
against H/3 produces a strongly orthogonal grid well suited for untangling the age-met allicity 
degeneracy. Our chosen Lick/IDS indices are summarized in Table [3l 



3.3. Emission-Line Correction 

The presence of H/3 emission superimposed on H/? absorption poses a significant prob- 
lem for our index measurements. The contamination of even weak H/3 emission lines in 
our galaxy spectra could significantly decrease the strengths of our measured H/5 absorption 
indices, thus systematically biasing our age distribution toward greater ages. Consequently, 
an emission correction was performed on our ho st-galaxy spectra. The basic procedure is 



outlined by lCaldwell. Rose. &: ConcannonI (120031 ). The method is based on two main asser- 



tions. First, although the underlying absorption-line spectrum at H/? is highly sensitive to 
age, the sensitivity to both age and metallicity is negligible at Ha. Second, the observed 
Hn/Fe index§| represents an effective indicator of the true underlying absorption-line spec- 
trum, implying that two galaxies with similar Hn/Fe indices will likewise possess similar Ha 
absorption strengths. 

Conservatively, we assumed that each galaxy in our sample was contaminated with 
Balmer emission from H II regions not included in the Vazdekis models. Therefore, we 
predicted the underlying Ha absorption spectrum for a given contaminated galaxy spectrum 
by matching that galaxy's Hn/Fe index with a Vazdekis model galaxy displaying a similar 



3Hn/Fe = <Hd7FeA4045 +H7/FeA4325 + H8/FeA3859>. 
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Hn/Fe index. The matching Vazdekis model SED was then treated as our model absorption 
spectrum. Next, we smoothed this spectrum to the appropriate resolution and velocity 
dispersion, normalized it to zero, and subtracted it from the host-galaxy spectrum. This left 
the Ha emission free of any contaminating absorption. The emission-line flux w as measured 
and converted into H/3 emission flux assuming a Case B Balmer decrement (jOsterbrock 



19891). 



We used a fiux-calibrated spectrum of the Orion Nebula taken with the same instrumen- 
tal setup to model the emission spectrum intended to be subtracted off the contaminated 
galaxy. The Orion spectrum was smoothed to the appropriate velocity dispersion and scaled 
so that the H/5 flux matched the calculated value. The continuum was removed from the 
modified Orion spectrum, and the result was subtracted from the contaminated host-galaxy 
spectrum. 



Results 



4.1. Diagnostic Grids 



In order to measure the H/5, FeA5270, and FeA4383 indices, we employed the FOR- 
TRAN77 code LECTOR provided by Alexandre Vazdekis. The code measures the indices of 
a one- dimensional input spectrum and calculates an error estimate b ased on photon statistics 
( ICardiel et al.lll998l : IVazdekis fc Arimotdll999l : ICenarro et al.ll200ll ) . With the line strengths 
we generated our index-index diagrams, in which each vertex represents a combination of 
line-index measurements for a single model of given age and metallicity. 

A qualitative comparison between several models and a representative host galaxy are 
given in Figure [T], showing a comparison of the spectrum for host galaxy NGC 4786 {Z ^ 
Zq) to model SEDs with similar metallicity {Z ^ Zq) but with varying population age. 
Quantitative age and metallicity estimates are found through a comparison with these model 
SEDs by way of an interpolation of the index-index diagnostic grids. Figure [2] shows the 
FeA5270-H/5 and FeA4383-H/5 grids. In each case the non-emission-corrected host-galaxy 
indices are plotted in the upper frames and the emission-corrected indices are plotted in the 
lower frames. Although for each galaxy metallicity and age interpolation we broadened the 
models according to its corresponding velocity dispersion, for the sake of presentation we 
show the whole sample plotted on a model grid broadened to a velocity dispersion of a = 
200 km s-^ 
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Fig. 1. — A comparison between host-galaxy NGC 4786 with an estimated metaUicity around 
solar and several Vazdekis model SEDs of varying ages. Each model has been broadened to 
the resolution and velocity dispersion of the host spectrum and labeled according to its age 
and H/9 index measurements. The series shows how the models change with age at a fixed 
metallicity and gives an alternate illustration of the technique used to realize our galaxy 
age/metallicity estimates. In this case, the age-sensitive index, H/3, for NGC 4786 most 
closely matches that measured for the SED at solar metallicity and a single population age 
of 7.08 Gyr. 



Computing ages and metallicities for our galaxy samples required the interpolation, 
and on occasion extrapolation, of an irregular gri d. This was ac c omplished us ing bivariant 
polynomial transformations of the following form (jWolbere]ll994j : ICardielll2003l ): 



N N-i 



agCj = y^^y^^ pjkjindxjy {indy 



k 

ij 1 
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N N-i 

Zi = qjk{indxiy{indyi)'', 

j k 

where N is the order of the polynomial and pjk and qjk are the coefficients that are solved 
for a set of nearest-neighbor points {mdx,indy) indexed by i. We performed a second-degree 
approximation, thus requiring us to compute at least 12 coefficients by solving two systems 
of 6 linear equations. An IDLq code was written which allowed us to interactively select the 
nearest 6-12 grid points (the number was dependent on the quality of the grid around the 
host galaxy's grid position) and the coefficients were computed by solving the above system 
of equations using the method of least squares. 

The upper and lower limits for age and metallicity were determined directly from the 
absorption-line index limits outputted by LECTOR. The extent of the index limits can be 
seen from the error bars in Figure [2l The upper and lower limits were treated as unique points 
for which age and metallicity estimates were made. In the case of large index errors, the grid 
points chosen by our code were different from those used to determine the age/metallicity of 
the galaxy. An assumption was made that the upper and lower H/3 index limits corresponded 
to the lower and upper limits, respectively, of our galaxy age estimate. Similarly, the upper 
and lower limit of the FeA5270 index corresponded to the upper and lower limits of the galaxy 
metallicity estimate. The final uncertainty was recorded as the magnitude of the difference 
between the age and metallicity of the limits and those of the galaxy. Furthermore, due to 
the irregular nature of the grid, a symmetric error in an index measurement translated into 
an asymmetric error in age and metallicity. 

The results from our analyses are shown in Table HI Column (1) gives the galaxy 
name. Columns (2)-(7) give the ages, metallicities, and corresponding errors determined via 
the FeA5270-H/5 diagram, while columns (8)-(13) give the same results ascertained via the 
FeA4383-H/3 diagram. 



"'Interactive Data Language: |http : / / www . itt vis . com/idl/idl7 . asp 
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4.2. Field Galaxy Sample 

Aspiring to test the uniqueness of SN la early-type host galaxies relative to random 
early-type galaxies, we compared the global properties of our SN la sample of elliptical host 
galaxies to those of a general sample of elliptical field galaxies from SDSS. The following 
steps were taken to generate a comparative sample of field galaxies. The SDSS Catalog 
Archive Server was queried for galaxy spectra within the redshift range < 2 < 0.15 that 
had a velocity dispersion measured by SDSS. The positions of these galaxies were then 
cross-referenced against galaxies in the NASA/IPAC Extragalactic Database (NED), and 
only those galaxies characterized as E/SO by NED were accepted. The galaxy spectra of 
these elliptical field galaxies were then obtained from the SDSS Data Archive Server. The 
spectra were emission-line corrected using a synthetic Orion spectrum created via the IRAF 
artificial data generation package, ARTDATA. Furthermore, rather than broadening the 
models to the specific velocity dispersion of each galaxy, we broadened the set of models to 
velocity dispersions in 30 km s~^ intervals from 100 km s~^ to 340 km s~^. The set of models 
used in the analysis for a given SDSS galaxy spectrum was the set broadened for a velocity 
dispersion that most closely matched that of the galaxy. Finally, to speed up the calculation, 
the grid points to perform the diagnostic interpolations were automatically chosen for the 
SDSS sample by an algorithm designed to mimic the grid-point selection criteria used for 
the interactive point selection in the host-galaxy analysis. 



4.3. Data Quality and Consistency Check 

The quality of our data is assessed via three checks, the results of which are shown 
in Figure [31 First, we compare the respective age and metallicity estimates obtained from 
our two index-index diagnostics. Figure [3^ shows the FeA4383 age estimate plotted against 
the age estimate of FeA5270. Given that both ages were determined using the same age- 
sensitive index, namely H/5, we would expect these data to be consistent with unity. As the 
top-left panel shows, we find moderate scatter about unity by virtue of our error estimates. 
Nevertheless, the data are consistent with unity, with a Xr of 0.88 for the fit. Figure[3]D shows 
the metallicity comparison between our two index-index diagnostics. In this case we see more 
significant dispersion; however, the data are still consistent with unity given the calculated 
uncertainty. Consequently, as a matter of preference, the galaxy ages and metaUicities used 
in the following analysis will be those determined from the H/3-FeA5270 diagnostic. 



Our second quality check is a comparison of o ur age and metallicity estimates to those 
of control galaxies studied by iTrager et al.l (I2OOOI ) (hereafter TOO). Our data are plotted 
against the data of TOO in Figure ^ and [3ll. TOO made age and metallicity measurements 
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2 2.4 2.8 3.2 2.4 3.6 4.8 6 
Fe5270 Fe4383 

Fig. 2. — Pre/post-emission-corrected diagnostic grids. Age-sensitive absorption-line in- 
dex equivalent width (A), H/9, is plotted against the metallicity-sensitive line indices EW, 
FeA5270 and FcA4383, before (top) and after (bottom) emission correction. The underlying 
grid reflect the Vazdekis model SEDs broadened to a velocity dispersion, o = 200.0 km s^^. 
Dashed lines (near vertical) are models of common metallicity, and dotted lines (near hori- 
zontal) are those of common age. The bold dashed line identifies models of solar metallicity, 
and the bold dotted hne marks models of age 7.08 Gyr. 
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Fig. 3. — Quality and consistency check. Panels (a) and (b) compare the age and metallicity 
estimates from our two diagnostic grids presented in Figure [2], respectively. In both cases we 
find minor variations that are nonetheless consistent with unity. Panels (b) an d (d) present 



acomparison showing reasonable agreement between our results and those of iTrager et al. 
J2000h . 
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through a central -Re/2 aperture and we are reporting the "preferred" non-solar abundance 
ratio (NSAR) model results (model 4). The age-sensitive line index used was H/5, while the 
two metallicity-sensitive indices were <Fe>|f| and Mg b. Figure [3t shows the age comparison 
for the two studies with the dotted line shows a model one-to-one relation. The plot shows 
fairly good agreement between the two studies. The xl of the one-to-one fit is 8.55, though 
it is improved if we throw out NGC 5813 at (8.69, 24.30), reducing xl to 2.72. 

The agreement is less encouraging for the metallicity comparison in Figure [311. Ignoring 
the few outlier points, the TOO galaxies seem to be systematically more metal rich than 
those measured in our study. This is likely due to the fact that we did not account for 
NSAR in our analysis. Studies have shown that the abundance ratios in early-type galaxies 
are often non-solar. In particular, the Mg/Fe abundan c e ratio has be e n shown to be larger 



in more luminous ea rly-type galaxies (lO'Connelll Il976l : iPeletierl Il989l : IWorthey et al.l Il992 



Vazdekis et al.lll997l ). Given that both Mg and Fe lines strongly factored into the metallicity 
measurements of TOO, we should not expect to see a straight 1:1 comparative ratio. Nev- 
ertheless, throwing out the points with the two greatest dispersions (NGC 5813 and NGC 
4489), we seem to be measuring the relative metallicities consistently with TOO, thereby 
enabling us to find potential trends between age/metallicity and SN la properties. 



We also compare our results to those published in lThomas et al.l (120051 ). In their study, 
they seek to set constraints on the epochs of early-type galaxy formation through an anal- 
ysis of 124 early-type galaxies in both high and low-density regions. We have 15 galaxies 
in common with their study and compare our H/5 index, age, and metallicity measurements 
to those from their analysis. Similar to the comparison with TOO, we see general agree- 
ment between our respective results. However, the agreement is far from one to one, with 
broad scatter greater than that seen in our TOO comparison. Once again, this not entirely 
unexpected given our lower S/N and the fact that we did not accounted for potential a en- 
hancement in our galaxies. Moreover, they were sampling a significantly smaller fraction of 
the galactic light w i th aii aperture radius of Re/10. Overall, our comparisons with both TOO 
and iThomas et al.l (120051 ) show reasonable agreement and confirm the ability of our analysis 
to measure acceptable relative population ages and metallicities for our SN la host-galaxy 
sample. 

Finally, we compare the B — K colors for the host galaxy and SDSS data to the galaxy 
age; see Figure Bl Th e expected average B — K color for a sample of elliptical galaxies is 
4.0 (jjarret et al.ll2003l ). and we can see that both our host galaxies and the SDSS galaxies 
meet this expectation. The figure also shows that both of the samples fall within the range 



5<Fe> = (i)(FeA5270 + FeA5335). 
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Fig. 4. — The SN la host-galaxy B — K color vs. age. We see that b oth of the samples 



fall •w ithin the range predicted by the evolutionary synthesis models of iBruzual fc Chariot 
( 120031 ). The center line is the predicted trend oi B — K color with age for a population of 
solar metallicity, the top line is the prediction for a metallicity of [M/H] ^ —0.68, and the 
bottom line is the prediction for a metallicity of [M/H] 0.40. 



predicted by the evolutionary synthesis models of iBruzual fc CharlotI (120031 ). The center 
line is the predicted trend of B — K color with age for a population of solar metallicity, the 
top line is the prediction for a metallicity of [M/H] ^ —0.68, and the bottom line is the 
prediction for a metallicity of [M/H] 0.40. 



4.4. Age and Metallicity 



The host-galaxy age and metallicity vs. SN la peak \^-band magnitude are plotted in 
Figure O The peak magnitudes were derived from measurements of the light-curve shape 
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parameter, A, fitted using MLCS2k2 by the authors of Ijha. Riess. fc Kirshnerl j2007l lR. In 



Figure [5^, the data show that the youngest populations tend to host the brightest SNe la, 
but with a wide spread. However, the eye has a tendency to gravitate toward large error bars, 
thereby placing an unwanted amount of attention to those points. Consequently, we have re- 
plotted the data in Figures [5]d and[5ll with point sizes that are inversely proportional to their 
uncertainty. In this way, the data with the smallest uncertainty in age will be represented 
by the largest points. When points with the largest uncertainty are de-emphasized, the 
trend between age and luminosity is very clear, but it is nevertheless difficult to distinguish 
whether the effect is a smooth transition with age or the result of two distinct populations. 



Gallagher et al.l (120051 ) found that only luminous SNe la occur in strongly star-forming 
hosts while E/SO galaxies show a wide range of SN la luminosities. This is confirmed in Figure 
[HI where we have plotted SN la peak magnitude vs. host-galaxy specific star formation, i.e., 
the SFR per stellar mas s. Th e SFRs were calculated from the Ha emission flux using 



the relation of iKennicuttI (119981 ). with the distance to each galaxy coming from the SN la 
luminosities. Solid circles represent the host galaxies from the current sample of early-type 
galaxies, while open circles are from the sample studied by G05. We eliminated galaxies 
from G05 if they were in our current sample since the Ha flux from G05 was uncorrected 
for underlying absorption. The points with arrows represent star formation upper limits 
in which Ha emission was buried within the noise. This trend, combined with Figure [Sjo, 
suggests that the age of the dominant population in E/SO galaxies determines the resulting 
SN la peak luminosity. A correlation bet ween SN la peak lu minosity and SFR has also 



been seen in the high-redshift SNLS data (jSuUivan et al.ll2006l ). but their type of search is 



biased against very low- luminosity events. Our "young" E/SO galaxies correspond to their 
"passive" hosts, and the extension to very old populations (with low specific SFR) confirms 
that population age is the major factor that determines the SN la peak luminosity. 

The plots of host metallicity vs. SN la peak luminosity (Figures and 01) also show 
a mild correlation. This may be due to the "age /metallicity degeneracy" shown in Figure 
[71 There we present measured galaxy age vs. metallicity for both the hosts and SDSS field 
galaxies, showing that the galaxies are not evenly distributed in the diagram but instead 
concentrated to the lower-left half. This is simply due to the evolution of the universe 
which started metal poor, meaning that there were few galaxies that are both old and metal 
rich and occupy the upper-right of the diagram. Thus, any correlation between age and 
a supernova property will result in some correlation with metallicity as well. Here, the 
relationship between host age and SN la luminosity is so clear that we claim it represents a 



6My(peak) = -19.504 + 0.736A + 0.182A2 + 5 log (Hq / 65 km s^^ Mpc^i) mag. 
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Fig. 5. — The SN la peak magnitudes (Hq = 72 km~^ s"-*^ Mpc~^) vs. (a) luminosity-weighted 
host-galaxy age and (c) metallicity. Panels (b) and (d) show the same data with the point 
size reflecting the average error along the abscissa. Point size is inversely proportional to 
the average error along the abscissa, with the exception of our young hosts in (b) where an 
upper limit is imposed on the point size. The data in panel (b) show a clear trend between 
SN la peak magnitude and host-galaxy age. The behavior could be evidence of either a 
smooth transition of SN magnitude with age or of two distinct populations of SNe la. Panel 
(d) shows a less convincing correlation likely arising from the degeneracy between age and 
metallicity. 
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Fig. 6. — The dependence of SN la peak magnitude on host-galaxy specific SFRho. Closed 
circles are the current sample of early-type hosts while open circles represent the high- 
SFR hosts of G05. The points marked by arrows indicate upper limits. On average, our 
specific SFR d i stribu tion is lower than that seen in the sample of high- 2; "passive" hosts of 
Sullivan et al.l (120061 ). However, aside from the expected decrease of SFR over cosmological 
time, it should be noted that our data only sample half of t he galactic light give n our 
extraction radius (-Re), and the SFRs for the "passive" galaxies of ISuUivan et al.l (l2006l ) were 
randomly assigned a rate of ~0.005 M0 yr~^. 
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Fig. 7. — An illustration of the age-metallicity degeneracy "problem" in the study of stellar 
populations. The plot shows the luminosity weighted age vs. metallicity results for our 
sample of SN la host galaxies and SDSS elliptical field galaxies. The plot reveals the trend 
stretching from young, metal rich to old, metal poor stellar populations. Histograms are also 
presented illustrating the age and metallicity distributions for the two respective samples. 



physical connection. We attribute mild effect in the metallicity relation to the lack of galaxies 
in the lower-right quadrant of Figure [Stl, which comes from age correlation combined with 
the fact that there are no old galaxies having high metal abundance in our sample. 

The distributions of age and metallicity for the two samples are also shown in Figure [TJ 
with the cumulative fraction plots shown in Figure [HI A Kolmogorov-Smirnov test reveals 
that there is a 63% chance that the host age distribution is drawn from the SDSS field galaxy 
age distribution, and a 44% chance that the host metallicity distribution is drawn from the 
field galaxy metallicity distribution. The results suggest that, given the size of the SN la 
sample, the abundance and age of early-type galaxies that host SNe la are similar on average 
to normal early-type field galaxies. Assuming this, the probability of a SN la going off in a 
given early-type galaxy does not strongly depend on the age or metallicity of the galaxy. 
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Fig. 8. — Cumulative fraction plots of the SN la host galaxy and SDSS galaxy age and 
metaUicity distributions. Our sample of SN la host galaxies is clearly similar in both age 
and metaUicity to the SDSS field galaxy sample. We find a 63% and 44% probability that the 
SDSS age and metaUicity distributions, respectively, represent acceptable parent populations 
from which the host metallicities were drawn. 
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4.5. Hubble Residuals 



The relation correcting SN la peak luminosity by the light-curve shape red uces the 
scatter in the low- 2; Hubble diagram to ~ 0.18 mag (jjha. Riess. fc Kirshnerl 120071 ) . If this 
residual scatter is purely random, then the evolution of age and metallicity in the universe 
will not induce a bias in cosmological measurements with SNe la. However, recent s tudies 
such as Timmes, Brown, & Truran (2003; hereafter TBT03) and iPodsiadlowski et al.l (120061 ) 
have predicted that SN calibration methods should be affected by a cosmic evolution in 
metallicity. 

In order to investigate this possible effect for SNe la hosted by early-type galaxies, 
we plot our measurements of age and metallicity against post-corrected SN la magnitude 
residual from the best fit to the Hubble diagram. The Hubble residual (HR) is defined as 

HR = /i^ — fJ'SN, 

where /i^ is the distance modulus determined from the host-galaxy redshiftQ, and fisN is the 
distance modulus determined via the MLCS2k2-corrected SN magnitude. HR is defined in 
the conventional manner such that overluminous SNe la have a negative Hubble residuajfj. 
The primary component to the uncertainty in HR comes from AfisN- The distance errors 
estimated from MLCS2k2 were found to be smaller than the scatter about the Hubble di- 



agram. The authors of iJha. Riess. &: Kirshnerl (120071 ) attributed this to an intrinsic scatter 
in SN la magnitudes and accordingly imposed a constant additive correction to these dis- 
tances. Given that we are probing for systematic changes in this intrinsic scatter with age 
and metallicity, we have not applied this correction to AfisN- 

The top-left panel of Figure [9] shows the Hubble residual vs. host-galaxy metallicity 
along with a least-squares fit to the data (solid line). The data and the linear fit show 
a correlation between host metallicity and deviation from the Hubble diagram such that 
supernovae in metal-rich galaxies are fainter than average and metal-poor galaxies host 
brighter than average events. The dashed line is a representation of the TBT03 relation 
(MtviVS. metallicity) showing that, at least qualitatively, the observed trend is consistent 
with that by expected by TBT03. 

To map the TBT03 prediction onto the plot, we converted the ^^ Ni yields from TBT03 
into V-band peak magnitude via the models of iHofiich et al.l (120021 ) as was done in Gal- 
lagher et al. (2005). To avoid an over-dependence on models, we also converted Ni masses 



^Imposed lower redshift cutoff at 2500 km s ^ to decrease uncertainty due to peculiar velocities. 
^Note that the opposite definition was adopted by Gallagher et al. ( 20051) . 
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Fig. 9. — SN la Hubble residual vs. luminosity-weighted metallicity (top left) and age 
(bottom left). A general trend is found suggesting that more metal rich populations produce 
intrinsically fainter SNe la. A least-squares fit finds to the data places a slope for the trend 
at 0.26. The dotted line represen ts the predicted trend based on the analytical models of 
Timmes. Brown, fc TruranI ( 120031 ). The results of our statistical tests are presented in the 
adjacent plots. We rule out a no-correlation result for the HR vs. metallicity data at the 
98% confidence level. A least-squares fit to the HR vs. age plot yields a nominal slope of 
—0.03 with negligible significance. 
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to peak magnitudes th rough an empirically determined relation using Ni mass estimated by 
Stritzinger et al.l (120061 ) and corresponding V-band luminosities given in lJha. Riess. &: Kirshner 
( 120071 ). Mapping TBT03 onto Figure using this conversion yields a qualitatively similar 
trend to both our data and the TBT03-Hoflich model. 

The top-right panel shows the results of the significance analysis for our linear fit. For 
each point ([M/H] ^3^-2)5 the Hubble residual was replaced by a randomly selected HR from 
a Gaussian distribution centered on the actual HR with a standard deviation equal to AHR. 
The same was done for each metallicity using the corresponding [M/H]^ and A[M/H]^. We 
performed 100,000 iterations, measured the best-fit slope (/3) to the new data each time, and 
generated the P distribution for the test. The results are seen in the top-right panel and 
show that a no-correlation result {(3 = 0) is ruled out at the 98% confidence level. A similar 
correlation between dust /extinction and Hubble residual is not seen for our sample, thereby 
removing the possibility of this correlation being a secondary dust effect. 

Assuming that the least-squares fit to the top-left panel of Figure [9] adequately repre- 
sents the trend of HR with metallicity, we would predict a 0.26 mag increase in intrinsic 
SN la magnitude per unit decline in progenitor (galaxy) metallicity. We have looked at sev- 
eral studies, two theoretical and one o bservational, tha t seek to ch aracterize c o smic chemical 
evolution. The theoretical studies by ICen fc Ostrikerl (119991 ) and ICen et al.l (120031 ) predict 
an approximate 0.05 dex and 0.25 dex drop in metallicity per un it redshift, respectively. The 
observational study conducted by iKewley fc KobulnickyI (|2007| ) find a preliminary result of 
0.15 dex/z. Taking the average of these three studies (0.15 dex/z) for the true decline in 
metallicity with redshift would translate into a 0.039 mag decrease in intrinsic SN la magni- 
tude per unit redshift. Assuming an ideal SN la data set out to z < 1.8, a minimization 
of three parameters (f^M, and hJ^ applied to a set of model SNe la subject to this 4% 
systematic error would induce an approximate 9% systematic error on the measurement of 
the equation-of-state parameter w = P/{pc^). 



Jha. Riess. fc Kirshnerl (120071 ) showed that the scatter of SN magnitudes about the 
Hubble diagram is lower in E/SO galaxies (~ 0.13 mag) than for sample as a whole (~ 0.18 
mag). Our sample of E/SO hosts show a scatter of approximately 0.14 mag. By applying 
a correction to SN magnitudes consistent with our least-squares fit, we reduce the scatter 
about the Hubble diagram to approximately 0.11 mag, illustrating the potential for such a 
correction to significantly improve SN la distance measurements. 

The bottom-left panel of Figure M shows the age of the parent galaxy vs. the Hubble 
residual of the SN la. The plot does not reveal a significant trend of HR with host-galaxy 



'Flat geometry assumed. 
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age. Indeed, the probability of a zero- correlation result is only ruled out at the 56% con- 
fidence level, suggesting that SN la magnitudes corrected for light-curve shape are likewise 
being corrected for any age bias. It should be noted that two galaxies (CGCG 016-058 and 
MCG-l-07-41-001) were removed from the age plot because their positions on their respective 
diagnostic grids required a difficult extrapolation that rendered their age estimates extremely 
uncertain. 



4.6. SN la Rate vs. Population Age 



Since the late 1970s, observations have shown that SNe la are more prevalent in star- 



pr 



forming late-type galaxies than in early- type galaxies (jOemler fc Tinsleyi Il979l ). This fact 
has been confirmed again and again with new studies showing that the SN la rate per 
unit mass is significantly higher in blue /late- type galaxies than in red /early-type galaxie s 
Jvan den Berghlll99ol : IPella Valle fc LiviolEggj : iMannucci et allbood : lleaman et aPboOsI ). 
The current explanation for these observations is that there are "prompt" and "delayed" 
( "tardy" ) SN la explosions. The prompt component is dependent on the rate of recent star 
formation, and the delayed component is dependent on the total number of low-mass stars. 
The combination of these two cornponents is belieyed to form the overall observed SN la 



rate (IScannapieco &: Bildstenll2005l : iMannucci et al.ll2005l : ISuUivan et al.ll2006l : iLeaman et al. 
2008h . 



Sullivan et al.l (120061 ) investigated the parameters shaping the overall SN la rate by 



observing the relative presence of SNe la in high- 2; host galaxies of differing mass and SFR. 
Our direct age measurement means that we can do a similar analysis at low redshifts by 
estimating the relative SN la rate as a function of host-galaxy age. Figures [TD] and [TT] 
summarize the details of our analysis on the relative SN la rate in elliptical host galaxies of 
various ages. 

The panel on the left of Figure [10] contains two unique histograms. The dotted-line 
histogram represents the number of host galaxies, or equivalently the number of SNe la, 
per age bin (identical to the distribution in Figure [7] with decreased bin size). The solid- 
line histogram is the expected number of stars per age bin. Each of the distributions has 
been normalized such that the bin encompassing 8 Gyr contains exactly one element. The 
number of stars is derived from the i^-band SDSS galaxy distribution, the K band being a 
relatively accurate tracer of stellar mass, i^'-band apparent magnitudes for the SDSS sample 
were compiled from NED and the distances to the galaxies in the sample were calculated 
from the SDSS redshift measurements, where care was taken to exclude those galaxies with 
Vr < 2500 km s~^. The i^-band absolute magnitudes, Mg, were derived from the distances 
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Fig. 10. — The relative SN la rate per stellar mass as a function of age in early-type galaxies. 
Stellar masses were approximated from each galaxy's i^-band luminosity. The left panel 
shows numbers of host galaxies, or SNe, per age bin along with the expected number of 
stars per age bin derived using the SDSS elliptical field galaxies. The right panel gives the 
expected relative SN la rate per age bin, showing a consistency with the "dual- component" 
model. 



and apparent magnitudes. Finally, the number of stars per galaxy was approximated by the 
expression 



The variable Mk,q is the corresponding absolute magnitude of the Sun in the K banco. 
The number of stars per age bin was then calculated by summing the stars in each galaxy 
in each age bin. 

The relative SN la rate per unit stellar mass is defined as the ratio of the number of SNe 
per age bin to the expected number of stars per age bin. The panel on the right of Figure [TUl 
shows the expected SN la rate per unit stellar mass as a function of age. The plot predicts 
a rate that is high for young galaxies, falls for intermediate-age galaxies, and increases and 
moderates for older populations. This result is consistent with the "dual-component model," 
given that the SN rate at low age could arise from a short delay time population of SNe la 
predicted to occur in galaxies with high SFR, and the rise in SN rate at late age could likely 



result from the predicted rise in SN rate with stellar mass (ISuUivan et al.ll2006l ). However 



'Mk^q = 3.28. 
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Fig. 11. — The relative SN la rate per stellar mass is replotted, showing the bimodal nature 
of the distribution. The dotted line represents the rate of white dwarf production while the 
dash-dot line is the accumulated number of white dwarfs as a function of time since the 
burst. 



we take the less general route and suggest the possibility that the rate at low population age 
is proportional to the rate of WD production following an initial burst of star formation. 
This function can be seen in Figure [TTl Similarly, the late-time SN rate can be traced by 
the cumulative number of WDs present in the galaxy (Figure [TT]) . This value is determined 
for an age, t = to, by integrating the WD production r ate from time t = t o t = to. These 



components of the SN la rate are related to those of ISuUivan et al.l (120061 ) . and yet they 
specifically target the population of stars believed to be the progenitors of SN la explosions. 

An important caveat to this analysis is that we are only sampling the dominant stellar 
population within a given host galaxy. It has been shown that a significant por tion (15 



30%) of elliptical and lent icular (SO) galaxies show evidence of recent star formation (lYi et al. 



2005 



Kaviraj et al.ll2007l ). which has the potential to compromise our conclusion that SNe la 



exploding in old stellar populations possess longer delay times. 
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5. Summary and Conclusions 

We have analyzed a sample of early-type galaxies that have hosted SNe la. Comparing 
the data with stellar population synthesis models, we have measured the global age and 
metallicity of the SN la elliptical host galaxies and a general sample of elliptical galaxies 
from the SDSS. Our results indicate that there is likely a significant correlation between 
age or metallicity and SN la absolute magnitude. Given the mounting evidence of a SN la 
rate dependence on specific SFR, we find it most likely that the SN la peak magnitudes are 
correlated with age, and the observed trend with metallicity is merely an artifact brought 
about by the evolutionary entanglement of age and metallicity. 

We find the global distributions of age and metallicity to be similar to those of our SDSS 
elliptical galaxy sample, suggesting that the presence, or absence, of a SN la in an elliptical 
galaxy is not dependent on either age or metallicity. Moreover, we detect a trend between 
early-type host-galaxy metallicity and the residuals from the Hubble diagram at the 98% 
confidence level. This trend is c onsistent with the pred ictions of TBT03 and with a trend 



observed for late-type galaxies by lGallagher et al.l (120051 ): it suggests that metal-rich galaxies 



produce underluminous SNe la, even after correcting for their light-curve shapes. Further- 
more, we conclude that the failure to apply a metallicity correction to SN la magnitudes 
could potentially introduce a 9% error into current and future measurements of w. 

Finally, we find the predicted SN la rate as a function of age for our sample of early-type 
ga laxies, and determin e it to be moderately consistent with the "dual- component model" 



of ISuUivan et al.l (120061 ). Also, we describe our preferred model governing the SN la rate in 
which the "prompt" component is proportional to the WD production rate and the "delayed" 
component scales with the cumulative WD population. 
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Table 1. Absorption- Line Galaxy Sample 



Galaxy SN My Morphology PA (°) Ap. Width (") L.C. ref.t 

(1) (2) (3) (4) (5) (6) (7) 



SN la host-galaxy Sample 



NGC 4374 


1991bg 


-17.88 


El 


90 


9.63 


1 


NGC 4526 


1994D 


-18.98 


SO 


110 


10.78 


1 


NGC 4493 


1994M 


-19.09 


E 


40 


4.03 


1 


CGCG 016-058 


1994T 


-18.67 


Sa 


-30 


3.01 


1 


NGC 2962 


1995D 


-19.42 


S0+ 





6.90 


1 


NGC 5061 


1996X 


-19.21 


Sa 


90 


4.17 


1 


NGC 5005 


1996ai 


-19.58 


Sbc 


-30 


12.66 


1 


NGC 5308 


1996bk 


-18.33 


SO- 


60 


10.78 


1 


NGC 5490 


1997cn 


-17.87 


E 





5.17 


1 


NGC 5440 


1998D 




Sa 


40 


6.19 




UGC 11149 


1998dx 


-19.14 


Sab 


-40 


6.75 


1 


NGC 7131 


1998CO 




SO 


-35 


4.45 


1 


NGC 6411 


1999da 


-17.81 


E 


70 


8.19 


1 


NGC 2841 


1999by 


-17.95 


Sab 


50 


5.88 


1 


NGC 6038 


1999CC 


-19.00 


Sc 


105 


8.19 


1 


NGC 2986 


1999gh 


-18.43 


E2;HII 


10 


6.90 


1 


UGC 11198 


2000dm 


-18.94 


Sab 


45 


7.06 


2 


CGCG 189-024 


2002G 


-19.347 


E 


10 


4.16 


2 


UGC 04322 


2002he 


-19.01 


E 


40 


5.32 


2 


CGCG 141-044 


2001bf 


-19.514 




100 


5.04 


2 


Anon 


2002aw 


-19.35 


Sb (f) 


80 


13.94 


2 


NGC 4786 


2002cf 


-17.82 


E-l-pec 


-20 


6.47 


2 


NGC 6702 


2002CS 




E 


55 


6.31 




MCG-F07-41-001 


2002do 


-18.48 


El 





20.85 


2 


NGC 7761 


2002ef 


-19.30 


SO 


90 


4.61 


2 


NGC 6986 


2002el 


-19.11 


SBO- 


5 


5.89 


2 


MCG-01-25-009 


2003D 


-18.31 


El 


10 


5.18 


2 


IC 0522 


2003Y 




SO 


-20 


7.34 




UGC 03787 


2003ch 




E-SO 





4.47 




NGC 4059 


2005bl 


-17.62 


E 


40 


4.61 


3 
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Table 1 — Continued 



Galaxy SN My Morphology PA (°) Ap. Width (") L.C. ref.t 
(1) (2) (3) (4) (5) (6) (7) 



Trager et al. (2000 ) Comparison Sample 



NGC 3608 • • • 


E2 


80 


6.19 


NGC 4472 • • • 


E2/S0 


-20 


12.80 


NGC 4478 • • • 


E2 


40 


6.17 


NGC 4489 • • • 


E 


110 


4.61 


NGC 4552 • • • 


E;HII 


-30 


6.31 


NGC 4649 • • • 


E2 


90 


13.09 


NGC 5638 • • • 


El 


-20 


8.04 


NGC 5813 • • • 


El-2 


-20 


6.17 


NGC 5846 • • • 


EO-1 


50 


11.35 


NGC 6127 • • • 


E 


90 


5.32 


NGC 6703 • • • 


SAO- 


90 


5.17 



Ul) Ijha. Riess. Kirshneil tom: ( 2) iGaneshalingam et al.1 (2008): 
FihDDenko et al.N200lh : iFiUDpenkol (|2005al) : (3) IXaubenbereer et al.1 (|2008l ): 
Gallagher et al.1 (|2005al ). 
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Table 2. Cross- Correlation Results 



Galaxy Vr (km s ^) 5vr (km s ^) a (km s ^) 



(1) 


(2) 


(3) 


(4) 


Anon 


7932.43 


19.66 


136.91 


CGCG 016-058 


10411.31 


17.79 


147.99 


CGCG 141-044 


4665.13 


23.01 


109.73 


CGCG 189-024 


10113.14 


21.01 


188.70 


IC 0522 


5092.26 


19.05 


170.76 


MCG-F07-41-001 


4584.86 


68.68 


344.72 


MCG-01-25-009 


6623.73 


21.58 


251.97 


NGC 2841 


633.49 


25.43 


244.72 


NGC 2962 


1967.99 


19.01 


186.61 


NGC 2986 


2310.81 


33.56 


283.67 


NGC 3608 


1237.80 


21.24 


209.54 


NGC 4059 


7195.99 


33.21 


210.88 


NGC 4374 


1020.48 


30.57 


316.56 


NGC 4472 


957.20 


30.95 


309.86 


NGC 4478 


1352.18 


17.67 


183.09 


NGC 4489 


948.10 


11.75 


106.89 


NGC 4493 


6957.26 


20.98 


214.39 


NGC 4526 


604.65 


24.90 


235.75 


NGC 4552 


316.64 


27.18 


282.82 


NGC 4649 


1087.03 


37.32 


344.45 


NGC 4786 


4611.97 


33.49 


288.86 


NGC 5005 


932.64 


23.29 


207.73 


NGC 5061 


2059.50 


22.15 


218.96 


NGC 5308 


2014.72 


22.08 


241.36 


NGC 5440 


3697.03 


27.35 


238.53 


NGC 5490 


4965.75 


37.73 


338.21 


NGC 5638 


1640.05 


20.15 


172.96 


NGC 5813 


1948.27 


28.57 


248.58 


NGC 5846 


1692.00 


27.71 


238.42 


NGC 6038 


9399.72 


32.83 


183.91 


NGC 6127 


4715.52 


29.99 


264.50 


NGC 6411 


3726.29 


23.52 


186.15 


NGC 6702 


4721.40 


32.71 


209.44 


NGC 6703 


2384.69 


20.62 


199.69 


NGC 6986 


8534.95 


25.52 


267.16 


NGC 7131 


5409.40 


36.77 


185.25 


NGC 7761 


7192.42 


19.90 


211.39 
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Table 2 — Continued 



Galaxy 


Vr (km s ^) 


6vr (km s ^) 


a (km s ^) 


(1) 


(2) 


(3) 


(4) 


UGC 11149 


16137.64 


92.35 


146.10 


UGC 11198 


4518.02 


19.91 


145.28 


UGC 03787 


8612.81 


57.54 


140.66 


UGC 04322 


7376.56 


23.98 


248.08 



Table 3. Lick/IDS Indices for the Study 





Index Bandpass 


Pseudocontinua 




Name 


(A) 


(A) 


Species Measured 


H/3 


4847.875-4876.625 


4827.875-4847.875, 4876.625-4891.625 


H, (Mg), (Cr), C 


Fe5270 


5245.650-5285.650 


5233.150-5248.150, 5285.650-5318.150 


Fe, C, (Mg), Ca 


Fe4383 


4369.125-4420.375 


4359.125-4370.375, 4442.875-4455.375 


Fe, C, (Si) 



Table 4. Ages and Metallicities 



Galaxy FeA5270 FcA4383 



[M/H] +<5 -<5 Age (Gyr) +5 -S [M/H] +S -<5 Age (Gyr) +5 -5 



(1) 


(2) 


(3) 


(4) 








\°} 


(9) 


(10) 


Cll^ 


('12") 




Anon 


0.081 




. ozo 


n 
u. 


ozo 


1.338 


0.194 


0.644 


0.158 




u. 


,^uo 




u. 


,oou 


1.338 


0.150 


0.644 


CGCG 016-058 


-0.374 


0. 


.278 


0. 


,266 


23.454 


13.895 


17.898 


-0.532 


0. 


,271 


Q 


.246 


26.083 


15.953 


22.769 


CGCG 141-044 


0.472 


0. 


.244 


0. 


,318 


1.273 


0.101 


0.101 


0.121 


0. 


,118 


0. 


.129 


1.256 


0.078 


0.061 


CGCG 189-024 


0.187 


0. 


.295 


0. 


,270 


2.611 


0.700 


0.694 


0.031 


0. 


,183 


0. 


.190 


2.891 


0.829 


2.027 


IC 0522 


-0.072 


0. 


.112 


0. 


,110 


8.200 


1.858 


2.849 


0.157 


0. 


,089 


0. 


,082 


6.583 


2.828 


3.035 


MCG+07-41-001 


-0.185 


0. 


.058 


0. 


072 


32.651 


5.932 


6.483 


-0.432 


0. 


,043 


Q 


,132 


36.121 


2.281 


1.886 


MCG-01-25-009 


-0.174 


0. 


.116 


0. 


107 


15.698 


5.170 


6.809 


0.041 


Q 


,100 


Q 


,081 


11.856 


3.307 


4.690 


NGC 2841 


0.037 


0. 


.069 


0. 


067 


12.558 


1.739 


2.606 


0.036 


0. 


,048 


Q 


,050 


12.578 


1.837 


3.073 


NGC 2962 


0.261 


0. 


.093 


0. 


,092 


5.180 


1.993 


2.185 


0.385 


0. 


,066 


0. 


,106 


4.210 


1.572 


1.818 


NGC 2986 


0.052 


0. 


.086 


0. 


,076 


9.652 


2.003 


2.108 


0.014 


0. 


,058 


0. 


,059 


9.667 


1.988 


1.998 


NGC 3608 


-0.013 


0. 


.067 


0. 


,069 


11.591 


2.072 


2.745 


0.011 


0. 


,050 


0. 


,052 


10.448 


1.235 


3.472 


NGC 4059 


-0.350 


0, 


.158 


0. 


,151 


7.567 


1.664 


3.247 


0.090 


0. 


,108 


0. 


,109 


5.514 


2.517 


2.934 


NGC 4374 


0.102 


0. 


.055 


0. 


052 


8.667 


1.195 


1.389 


0.081 


0. 


,045 


Q 


,048 


8.851 


1.369 


1.700 


NGC 4472 


0.179 


0. 


.044 


0. 


,045 


8.752 


1.273 


1.466 


0.154 


0. 


,041 


0. 


,038 


9.360 


1.325 


1.346 


NGC 4478 


-0.088 


0. 


.079 


0. 


,077 


9.027 


1.641 


1.964 


-0.031 


0. 


,051 


0. 


.057 


8.838 


1.882 


1.812 


NGC 4489 


0.428 


0. 


.186 


0. 


200 


2.539 


0.719 


0.753 


0.167 


0. 


,099 


0. 


.104 


3.350 


0.808 


2.099 


NGC 4493 


0.014 


0. 


.221 


0. 


,224 


7.165 


3.964 


5.585 


-0.077 


0. 


,161 


0. 


,160 


7.561 


4.088 


6.693 


NGC 4526 


0.400 


0. 


.061 


0. 


062 


1.669 


0.215 


0.197 


0.488 


0. 


,036 


Q 


,031 


1.122 


0.279 


0.312 


NGC 4552 


-0.040 


0. 


.024 


0. 


063 


13.481 


1.762 


2.532 


0.056 


0. 


,047 


Q 


,041 


11.834 


1.136 


2.511 


NGC 4649 


-0.076 


u. 


.Uoi 


u. 




21.565 


2.507 


2.533 


0.089 


u. 




u. 


,U4o 


16.931 


2.295 


2.219 


NGC 4786 


0.000 


0, 


.099 


0, 


,097 


7.100 


1.684 


1.740 


0.092 


0, 


,077 


0, 


,067 


6.215 


2.690 


2.219 


NGC 5005 


































NGC 5061 


0.359 


0, 


.035 


0. 


,056 


1.923 


0.185 


0.054 


0.427 


0, 


,026 


0, 


.026 


1.429 


0.271 


0.202 


NGC 5308 


-0.177 


0, 


.063 


0, 


,063 


11.087 


1.923 


1.809 


0.025 


0, 


,054 


0, 


,053 


9.578 


1.758 


1.934 


NGC 5440 


0.057 


0, 


,107 


0, 


102 


7.415 


1.957 


2.536 


0.248 


0, 


,076 


0, 


,080 


5.453 


2.151 


2.645 


NGC 5490 


0.012 


0, 


,124 


0, 


120 


9.840 


2.867 


3.296 


0.155 


0, 


,099 


0, 


,097 


8.113 


3.096 


3.338 


NGC 5638 


-0.105 


0, 


,079 


0, 


,081 


13.743 


3.146 


3.677 


0.175 


0, 


,054 


0, 


,060 


9.903 


2.322 


2.165 


NGC 5813 


0.027 


0, 


,086 


0, 


,074 


8.692 


1.891 


1.996 


0.369 


0, 


,068 


0, 


,050 


4.921 


0.737 


2.135 


NGC 5846 


-0.147 


0, 


,099 


0, 


,047 


14.421 


2.525 


4.724 


-0.038 


0, 


,045 





.053 


13.536 


2.046 


3.490 


NGC 6038 


-0.302 


0, 


,194 


0, 


,191 


12.286 


5.161 


10.610 


-0.482 


0, 


,186 


0, 


.167 


13.933 


6.198 


8.528 


NGC 6127 


-0.178 


0, 


,121 


0, 


,066 


14.917 


3.711 


6.099 


-0.122 


0, 


,097 


0, 


.050 


13.975 


3.287 


5.596 


NGC 6411 


-0.108 


0, 


,092 


0, 


,091 


14.833 


3.788 


5.096 


-0.442 


0, 


,067 





.067 


19.433 


4.681 


5.155 


NGC 6702 


0.442 


0, 


,121 


0, 


,126 


1.965 


0.379 


0.530 


0.157 


0, 


,082 


0, 


,068 


2.845 


0.478 


0.548 


NGC 6703 


0.078 


0, 


,069 


0. 


,063 


8.186 


1.451 


1.695 


0.213 


0, 


,052 


0, 


.052 


6.829 


1.805 


1.727 


NGC 6986 


0.602 


0, 


,264 


0, 


,320 


1.289 


0.111 


0.067 


0.744 


0, 


,225 


0, 


.122 


1.020 


2.512 


0.179 



Table 4 — Continued 



Galaxy FeA5270 FeA4383 





[M/H] 


+S 


-S 


Age (Gyr) 


+S 


-S 


[M/H] 


+S 


-S 


Age (Gyr) 


+S 


-S 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


NGC 7131 


-0.281 


0.134 


0.130 


23.653 


7.315 


8.609 


-0.327 


0.112 


0.107 


24.228 


7.778 


9.103 


NGC 7761 


0.164 


0.186 


0.192 


2.907 


0.613 


0.803 


0.268 


0.128 


0.085 


2.670 


0.652 


0.864 


UGC 11149 














-0.480 


0.251 


0.245 


9.768 


5.339 


12.562 


UGC 11198 


-0.238 


0.108 


0.112 


18.449 


6.004 


5.519 


-0.204 


0.085 


0.083 


17.861 


5.797 


5.749 


UGC 03787 


0.446 


0.291 


0.419 


1.284 


0.076 


0.082 


0.422 


0.186 


1.729 


1.172 


0.235 


0.190 


UGC 04322 


-0.155 


0.155 


0.137 


11.844 


3.818 


6.540 


-0.116 


0.108 


0.102 


11.527 


3.927 


6.212 



